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Abstract 



o 
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I ' Theory of melting of the vortex lattice in type II superconductors in the 

Q^' framework of Ginzburg - Landau approach is presented. The melting line 

C/3 



location is determined and magnetization and specific heat jumps along it are 
calculated . The magnetization of liquid is larger than that of solid by 1.8% 
irrespective of the melting temperature, while the specific heat jump is about 
6% and decreases slowly with temperature. The magnetization curves agrees 
'^ ■ with experimental results on YBCO and Monte Carlo simulations. 
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I. INTRODUCTION AND THE MAIN IDEA 

Abrikosov vortices created by magnetic field in type II superconductors strongly interact 
with each other creating highly correlated configurations like the vortex lattice. In high Tc 
cuprates thermal fiuctuations at relatively large temperatures are strong enough to melt 
the lattice. Several recent remarkable experiments discovered that the vortex lattice 
melting in high T^ superconductors is first order with magnetization jumps fl]] and spikes 
in specific heat [0. Magnetization and entropy jumps were measured using great variety of 
techniques: local Hall probes ^, SQUID |0^, torque magnetometry f^J^ and integrating 
the specific heat spike [^. It was found that in addition to the spike there is also a jump in 
specific heat which was measured as well p|,|3l- These precise measurements pose a question 
of accurate quantitative theoretical description of thermal fiuctuations in vortex matter. 
The Ginzburg - Landau (GL) approach is appropriate to describe thermal fiuctuations 
near T^. The GL model is however highly nontrivial even within the lowest Landau level 
(LLL) approximation valid at high fields. In this simplified model the only parameter is 
the dimensionless scaled temperature a^ ~ (T — T^f (H)) / {T HY^^ (defined more precisely 
in eq.( p.2a| ) below). Over last twenty years great variety of theoretical methods were 
applied to study this model. Brezin, Nelson and Thiaville [Q applied the renormalization 
group (RG) on the one loop level. No fixed points of the (functional) RG equations were 
found and they concluded therefore that the transition from liquid to the solid is first order 
0. This approach however does not provide a quantitative theory of melting. 
Two perturbative approaches were developed and greatly improved recently to describe the 
solid phase and the liquid phase respectively. The perturbative approach on the liquid side 



was pioneered by Thouless and Ruggeri |T0[ . They developed a perturbative expansion 
around a homogeneous (liquid) state in which all the "bubble" diagrams are resummed. 
Unfortunately they found that the series are asymptotic and although first few terms 
provide accurate results at very high temperatures, the series become inapplicable for ot 
less than —2 which is quite far above the melting line (believed to be located around 
ot ~ —10). We recently obtained the optimized gaussian series [|ll| which are convergent 
rather than asymptotic with radius of convergence of a^ = —5 (a^ = —5 is still above the 
melting point). 

On the solid side, Eilenberger and Maki and Takayama ||T^ calculated the fiuctuations 
spectrum around Abrikosov's mean field solution. They noticed that the vortex lattice 
phonon modes are softer than that of the acoustic phonons in atomic crystals and this leads 
to infrared divergences in certain quantities. This was initially interpreted as destruction 
of the vortex solid by thermal fiuctuations and the perturbation theory was abandoned. 
However the divergence looks suspiciously similar to "spurious" IR divergences in the 
critical phenomena theory and recently it was shown that all these IR divergences cancel in 
physical quantities |TB|. The series therefore are reliable, and were extended to two loops, 
so that the LLL GL theory on the solid side is now precise enough even at melting point. 
However on the liquid side one needs a theory in the region —10 < a^ < —5. Moreover this 
theory should be very precise since free energies of solid and liquid differ only by few 
percents near melting. This requires a better qualitative understanding of the metastable 
phases of the theory. It is clear that the overheated solid becomes unstable at some finite 
temperature. It not clear however whether overcooled liquid becomes unstable at some 



finite temperature (like water) or exists all the way down to T = as a metastable state. 
The gaussian (Hartree - Fock) variational calculation, although perhaps of a limited 
precision, is usually a very good guide as far as qualitative features of the phase diagram 
are concerned. Such a calculation in the liquid was performed long ago |]10|, while more 
complicated one sampling also inhomogeneous states (vortex lattice) was obtained recently 
||TT| . The gaussian results are as follows. The free energy of the solid state is lower than 
that of the liquid for temperatures lower than melting temperature. The solid state is 
therefore the stable one below it, becomes metastable at somewhat higher temperatures 
and is destabilized at a^ = —5. The liquid state becomes metastable below the melting 
temperature, but unlike the solid, does not loose metastability all the way down to T = 0. 
The excitation energy of the overcooled liquid approaches zero as a power e ~ 1/a^. 
Meantime in different area of physics similar qualitative results were obtained. It was 
shown by variety of analytical and numerical methods that liquid (gas) phase of the 
classical one component Coulomb plasma exists as a metastable state down to T = with 
energy gradually approaching that of the Madelung solid and excitation energy 
diminishing |T^. It seems plausible to speculate that the same would happen in any 
system of particles interacting via long range repulsive forces. In fact the vortices in the 
London approximation become a sort of repelling particles with the force even more long 
range than Coulombic. This was an additional strong motivation to consider the above 
scenario in vortex matter. Below we provide both theoretical and phenomenological 
evidence that the above scenario is the correct one. 

Assuming absence of singularities on the liquid branch allows to develop an essentially 
precise theory of the LLL GL model in vortex liquid (even including overcooled liquid) 
using the Borel - Fade (BF) ||15| method at any temperature. First we clarify several issues 



which prevented the use and acceptance of the BF method in the past and then combine it 
with the recently developed LLL theory of solids to calculate the melting line and 
magnetization and specific heat jumps. Very early Ruggeri and Thouless |]TD| tried to use 



BF to calculate the specific heat without much success because their series were too short 
jlGl . Attempts to use BF for the calculation of melting also ran into problems. Hikami, 



Fujita and Larkin |T^ tried to find the melting point by comparing the BF energy with the 
one loop solid energy and obtained aj- = —7. However their one loop solid energy was 
incorrect and in any case it was not precise enough. 
The LLL GL model was also studied numerically in both 3D [l^ and 2D ||l9l and by a 



variety of nonperturbative analytical methods. Among them the density functional 
1/A^ [^, dislocation theory of melting |2^ and others |2^ . 



As we show in this paper, the BF energy combined with the correct two loop solid energy 
computed recently gives scaled melting temperature a™ = —9.5 and in addition predicts 
other characteristics of the model. The melting line location is determined and 
magnetization and specific heat jumps are calculated . The magnetization of liquid is 
larger than that of solid by 1.8% irrespective of the melting temperature, while the specific 
heat jump is about 6% and decreases slowly with temperature. 

In addition to theory of melting we calculated magnetization and specific heat curves. The 
magnetization curves agree quite well with Monte Carlo simulations of the LLL GL p 



and almost perfectly for specific heat in 2D by Kato and Nagaosa in ref. |19|. However to 



describe experimental results like YBCO at lower fields and lower temperature, higher 



Landau levels (HLL) corrections are required. Experimentally it was claimed that one can 



establish the LLL scaling for fields above 3T 1^^. A glance at the data however shows that 



above T^ the scaling for magnetization curves is generally very bad: the magnetization due 
to the LLL contribution is much larger that the experimental one above T^. Therefore we 
calculated the leading correction by the HLLs and then compare with experiments. 
The paper is organized as follows. The models are defined in section II and the melting 
theory in the LLL is described in section III. In section IV, the HLL correction is discussed 
and the magnetization curves are compared with experiments. We conclude in section V. 

II. THE MODEL AND BASIC ASSUMPTIONS 

A. The model 

Thermal fiuctuations of 3D materials with not very strong asymmetry along the z axis are 
effectively described by the following Ginzburg-Landau free energy: 

where A = [By^ 0) describes magnetic field (considered constant and nonfiuctuating, see 

below) in Landau gauge and covariant derivative is defined by 

D = V — 2^A,$o = ^(e* = 2e). Statistical physics is described by the statistical sum: 

D^D^exp ("-0 (2) 

Our aim is to quantitatively describe the effects of thermal fluctuations of high T^ cuprates 
of the YBCO type on the few percent precision level (such a precision is required for the 
theory of melting since energies of liquid and solid near melting differ only by a few 
percent). 



B. Assumptions 

The use of the above GL energy makes several physical assumptions. They are listed below. 

(1) Continuum model 

We use anisotropic GL model despite the layered structure of the high T^ cuprates for 
which the Lawrence - Doniah model is more appropriate model. Effects of layered 
structure are dominant in BSCCO or Tl compounds ( 7 = \frri^m^\j > 1000) and 
noticeable for cuprates with anisotropy of order 7 = 50 like LaBaCuO or Hgl223. The 
requirement, that the GL can be effectively used, therefore limits us to optimally doped 
YBCOj-s (or slightly overdoped or underdoped) for which the anisotropy parameter is not 
very large 7 = 4 — 8), DyBCO or possibly Hgl221 which has a slightly larger anisotropy. 

(2) Range of validity of using the mesoscopic (GL) approach 

The GL approach generally is an effective mesoscopic approach applicable when we can 
neglect higher order terms generated when one "integrates out" microscopic degrees of 



freedom typically but not always near second order phase transitions. The leading higher 
dimensional terms we neglect (as "irrelevant") are \il)\^ and higher (four) derivative terms. 
After that the model practically becomes rotationally symmetric in the ah plane. We can 
rescale it to ma — mi, = m,ab- For several physical questions this assumption is not valid 
because irrelevant terms neglected might become " dangerous" . For example the question 
of the structural phase transition into the square lattice is clearly of this type p5[. It is 
known that even assuming ma/m,b = 1 in low temperature vortex lattices in YBCO 
rotational symmetry is broken down to the fourfold symmetry by the four derivative terms. 
However there is no significant correction to the magnetization from the those higher 
dimension terms. 

(3) Expansion of parameters around T^ 

Generally parameters of the GL model of eq. (|I]) are complicated functions of temperature 
which is determined by the details of the microscopic theory. We expand the coefficient 
a{T) near Tc : 

a{T) = Tja(l - t) - a'(l - tf + ...], (3) 

where t = T/T^. The second and higher terms in the expansion are omitted and therefore 
when temperature deviates significantly from T^ one cannot expect the model to have a 
good precision. 

(4) Constant nonfluctuating magnetic field 

For strongly type II superconductors like the high T^ cuprates not very far from Hc2{T) 
(this easily covers the range of interest in this paper, for the detailed discussion of the 
range of applicability beyond it see ref. p^) magnetic field is homogeneous to a high 
degree due to superposition from many vortices. Inhomogeneity is of order I/k^ ~ 10~^. 
Since the main subject of this study is the thermal fiuctuations effects of the order 
parameter field, one might ask whether thermal fiuctuations of the electromagnetic field 
should be also taken into account. Halperin, Lubensky and Ma considered this question 
long ago 1^. The conclusion was that they are completely negligible for very large k. 



Upon discovery of the high Tc cuprates, the issue was reconsidered |^ and the same result 
was obtained to a very high precision. Therefore here magnetic field is treated both as 
constant and nonfiuctuating {B = H) and the last term in eq.(P can be omitted (to 
precision of order 1/k^). However when we calculate the magnetization, M = [B — H)/A7r 
which is of order I/k^, high order correction must be considered. 
(5) Disorder. 

Point - like disorder is always present in YBCO. For example when the melting line 
becomes the second order transition and magnetization becomes irreversible the systems 
show clear disorder effects. However in some samples the disorder effects are minor. In the 
maximally oxidized sample H the second order transition is not seen even at highest 



available fields (30T). In the optimally doped sample the same is true up to 12T |2^. The 
situation is believed to be different in other materials like DyBa2CusOj [fTlpOf or twinned 
YBCO. We will address the disorder problem in our future publication. 



C. Landau level modes in the quasimomentum basis 

Assuming that all the requirements are met we now divide the fluctuations into the LLL 
and HLL modes. Throughout most of the paper will use the coherence length 
^ = \/f? I {27nabOiTc) as a unit of length, Tc as unit of temperature, T = tTc, and 
— "dT (^ ~ 2ff2 ^^ ^ ^^^^ °^ magnetic field, B = bHc2- After rescaling eq.(pby 

X ^ ^x,y ^^y,z^ ^,ip'^ ^ ^ip'^ (7 = ^/mjm^b) one obtains: 
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It is convenient to expand the order parameter field in a complete basis of noninteracting 
theory: Landau levels. In the hexagonal lattice phase the most convenient basis is the 
quasimomentum basis: 
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Here (p'^{x) the n Landau level with quasi- momentum k: 
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Even in the liquid state which is more symmetric than the hexagonal lattice we find it 
convenient to use this basis. 



III. MELTING LINE, MAGNETIZATION AND SPECIFIC HEAT IN LLL 

APPROXIMATION 



The LLL contributes significantly the thermal fluctuation which in particular leads to the 
melting (for example, ref.( [§|), and Pierson and Vails in ref. ( p^) et at). In this section, 
we will essentially solve the LLL GL model. 



A. The LLL scaling 

If the magnetic field is quite high, we can keep only the n = LLL modes (in eq.(^, we 
use the LLL condition jD^/'p = 6|'?/'p ): 



■' T u 



d X 



l|9.^p-a.H^ + i|# + ^li^^ 



(9) 



Within the LLL approximation the problem simplifies (no gradient term in directions 
perpendicular to the field) and possesses an LLL scaling |^. After additional rescaling 

/ \-l/3 / \2/3 

X -^ x/vb, y -^ y/Vb, z ^ z i -—y= j , ^^ — > ( -^ j ip"^, the dimensionless free energy 



becomes: 
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and it can be approximated to order 4 
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where the scaled temperature 



Free energy density in the newly scaled model is: 
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is a function of a^ only to order ^. 



(10) 
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B. Free energy of liquid and solid 

Now we specify the solution of the LLL GL model. The liquid LLL (scaled) free energy is 
written as 



The function g can be expanded as 



(14) 



,.x 



(15) 



where the high temperature small parameter x = \e ^' ^ is defined as a solution of the 
Gaussian gap equation 
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(16) 



for the excitation energy e. The coefficients c„ can be found in |T^. We will denote by 
gk {x) the [/c, k — 1] BP transform [|l5l of g{x) (other BP approximants clearly violate the 
correct low temperature asymptotics). The BP transform is defined as 
Jq°° g',^ (xt) exp (— t) dt where g[ is the [k, k — 1] Pade transform of "^^=1 ^^• 
For k = 4, and k = 5, the liquid energy converges to required precision (0.1%). In what 
follows we will use g^. 
The solid energy to two loops is [0,R2]: 
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(17) 



where (3a = 1.1596. On Fig.l we plot the energies of solid and liquid. They are very close 
(see the difference on inset). The liquid energy completely agrees with the optimized 
gaussian expansion results [^ till its radius of convergence at ot = — 5. The energy of the 
overcooled liquid at low temperatures approaches that of the vortex solid. We obtained 



similar results in 2D and found that they also agree with Monte Carlo simulations [^ and 
optimized gaussian expansion W^ (see a brief discussion in the next section) 



C. Melting line. Comparison with experiments, Monte Carlo simulations and 

Lindemann criterion 

Comparing solid and liquid energy (inset in Fig.l), we find that a^ = —9.5. This is in 
accord with experimental results. As an example on Fig. 2 we present fitting of the melting 
line of fully oxidized YBa2Cu30Y [| which gives T^ = 88.2, H^2 = 175.9, Gi = 7.0 10"^ 
Melting lines of optimally doped untwinned 0,^,0 YBa2CusOi^s and DyBa2Cu307 
are also ffited extremely well. The results of ffiting are given in Table 1. 



Table 1 

Parameters of high Tc superconductors deduced from the melting line 



material 


T 


Hc2 


Gi 


K 


7 


reference 


YBCOj-s 


93.07 


167.53 


1.910-^ 


48.5 


7.76 


M 


YBCOj 


88.16 


175.9 


7.010-5 


50 


4 


i 


DyBCOe.7 


90.14 


163 


3.210-5 


33.77 


5.3 


i 



The available 3D Monte Carlo simulations ||18[ unfortunately are not precise enough to 
provide an accurate melting point since the LLL scaling is violated and one gets values of 
a^ = —14.5, —13.2, —10.9 at magnetic fields 1, 2, 5T respectively. We found also that the 
theoretical magnetization calculated by using parameters given by ref. [|I^ is in a very good 
agreement with the Monte Carlo simulation result of ref. [|18| . However the determination 
of melting temperature needs higher precision, and the sample size (~100 vortices) used in 



ref. [^ may be not large enough to give an accurate determination of the melting 
temperature (due to boundary effects, LLL scaling will be violated too). The situation in 
2D is better since the sample size is much larger. We performed similar calculation for the 



2D LLL GL liquid free energy, combined it with the earlier solid energy calculation [13,32 



/.o. = -J^ + 21og^-lM_2.92. (18) 

and find that the melting point a™ = —13.2. It is in good agreement with MC simulations 

El- 

Phenomenologically melting line can be located using Lindemann criterion or its more 
refined version using Debye - Waller factor. The more refined definition is required since 
vortices are not pointlike. It was found numerically for Yukawa gas [Q that the Debye - 



Waller factor e^^^ (ratio of the structure function at the second Bragg peak at melting to 
its value at T = 0) is about 60%. Using methods of ||3^, one obtains for the 3D LLL GL 
model 



e-^"^ = 0.59. (19) 



D. Magnetization jump at melting 

The scaled magnetization is defined by: 

m M = --^feff (ot) (20) 

At the melting point a^ = —9.5 the magnetization jump AM divided by the 
magnetization at the melting on the solid side is 

^ = ^ = .018 (21) 

Ms nis 

This is compared on the inset of Fig. 2 with experimental results of fully oxidized 
YBa2Cu^0j (rhombs) and optimally doped untwinned YBa2Cu^Oj^5 [Q (stars). The 
agreement is quite good. 



E. Specific heat jump at melting 

In addition to the delta function like spike at melting for specific heat experiments shows 
also a specific heat jump. The theory allows to quantitatively estimate it. 
The specific heat of the vortex lattice is C = —T-^F(,ff. The scaled specific heat is 
defined as 



c 160 A f huj y^^ 

V4vrV2/ 

3t2 ^ ^ V47rv^y 

.|i(2-26 + t)V"(aT) 



+ '-^ih-l-t)(J^T fM (22) 



Hl,T 



where Cmf is the mean field specific heat of solid, ^^^jg ^2 ■ 
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Thus the specific heat jump is: 

Ac = 0.0075 ( '^~^l^^^ ) -0.20Gi'/^{b-l-t)f^] (23) 



Using the parameters obtained fitting the melting hne, we compare eq.(P3D with the 
experimental result of ref. M in the inset of Fig. 2. 



IV. HIGHER LANDAU LEVELS CONTRIBUTIONS 

A. Where is the LLL approximation really valid? 

Contributions of HLL are important phenomenologically in two regions of the phase 
diagram. The first is at temperatures above the mean field Tc inside the liquid phase. The 

second is far below the melting point deep inside the solid phase. 

Naively when "distance from the mean field transition line" is smaller than the "inter 
Landau level gap" , 1 — t — b < 2b one expects that higher Landau harmonics can be 
neglected. More careful examination shows that a weaker condition 1 — t — b < 126 should 
be used for a validity test of the LLL approximation [^ to calculate the mean field 



contributions in vortex solid. Additional factor 6 comes from the hexagonal symmetry of 
the lattice since contributions of higher Landau levels (HLL) , first to fifth HLL do not 
appear in perturbative calculation. How what about contributions beyond mean field? The 



question has been studied by Lawrie ||35[ in the framework of the Hartree - Fock (gaussian) 
approximation in the vortex liquid. The result was that the region of validity is limited. In 
this section we will incorporate the leading HLL correction using gaussian approximation 
and then compare the theoretical results with experimental magnetization curves. 



B. Gaussian Approximation in the liquid phase 

The free energy density is: 

Feff = -^^InJDijDi^expi-f) (24) 

where /, V is the scaled full GL model eq.(||) and the scaled volume. In gaussian 
approximation, / is divided into an optimized quadratic part K, and a "small" part v. 
Then K is chosen in such a way that the gaussian energy is minimal. The gaussian energy 
is a rigorous lower bound on energy. Due to translational symmetry of the vortex liquid an 
arbitrary symmetric quadratic part K has only one variational parameter e, 

^ = ;^ (^ (iD^r - &i^r) + l\d.^\' + ^i^r) (25) 

The small perturbation is therefore: 
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{-a,-e)\^\' + -\^\' 
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The gaussian energy consists of two parts. The first is the Trlog term: 



27rfi;2 



log 



Vtpexp^-K) 



uH% 1 



'c2 



27rK2 



^ hy y/nb + e, 

7277 ^ 
* n=0 



(27) 



The second is proportional to the expectation value of f in a solvable model defined by K 

2" 






[v) 



coH^ 



c2 



2ttk'^ 



-ah - £) 






2V27f^^Vnb + 6 



+ UJ 



2v/2' 



vr 






n=0 



-\/n6 + £ 



(28) 



Both are divergent in the ultraviolet. Introducing a UV momentum cutoff which effectively 



limits the number of Landau levels to Nj = j — 1, the Trlog term diverges as: 



^ by Vnb + e = ^ 

* n=0 



2^3/2 + (. - ^ ] K-l- 



+ ftrlog{£,b) 



(29) 



with the last term finite. The "bubble" integral diverges logarithmically: 



(30) 



'^H^c2 ' 



Substituting eq. (|29|) into the gaussian energy one obtains: (in units of ■^:^ 
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A'/' + ^(-l-Ai/2)2 + (^_a, - -L_Ay2 _ a,dJ,rio,{e, b) + 2u^K"'dJtrXog{e, b) (31) 



7277 3 

-edeftr\og{e, b) + uj{dJtriog{e, b)f + ftr\og{e, b). 

Tc will be renormalized. We define a^ = a'f^ + 2uj^-A^^^, the above energy becomes (in 



V27I 



units of TT-^' 



1 .^A3/2_ ^^^^1/2)2 



—ay 



1 



7277 3 '7277 ' V 2; 7277 

-edeftr\og{e, b) + Uj{deftr\og{e, b)f + ftr\og{e, b). 



A'/^-aldJtriogie,b) 



(32) 



Thus the temperature and vacuum energy will be renormalized. The first three terms are 
divergent and linear in temperature, they will not contribute to any physical quantities, like 
magnetization and specific heat. Minimizing the energy eq. (|32|) , we get the gap equation: 

e = -al - 2udJtT\og{s, b) (33) 

Superscript "r" will be dropped later on. Function ftriog{£,b) can be written 

ftriog{e,b) = ^b'/'g[^ (34) 
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where 



°^ r 9 19 1 

9{x) = Y^ Vn + X - -{x + n + -)^ + -{x + n - - 



n=0 '- 



l^^-¥ ^^^^ 



For the LLL model in the gaussian approximation, g (x) = \px. For Prange hmit |3^ 
Gi -^ 0, the free energy will be '^^h^''^g {-f). 



C. Magnetization 



When K, is quite large and magnetization can be approximated by 

M = -^F,ff{B). (36) 

The HLL correction will be calculated as follows. We numerically solve the gap equation 
(0) , and Feff{B) can be obtained. Then eq.(^) is used to calculate the magnetization of 
the full GL model in gaussian approximation. The HLL correction is thus the 
magnetization of the full GL model in gaussian approximation minus the magnetization of 
the LLL contribution in gaussian approximation. We compare the experiments using 



following approximation. The LLL contribution (for example, see ref.( |jTT|)) is 



MLLL = ^,-m{aT), (37) 

where m (ot) is defined in eq.(^) is taken exactly (calculated by the BP method to 
highest order), while the corrections due to HLL are taken in gaussian approximation. The 
comparison of the theoretical predictions with the experiments for fully oxidized 
YBa2CuzOj Q, is shown on Fig. 3. We use the experimental asymmetry value 7 = 4 and 
values of T^, Hc2 and Gi from the fitting of the melting curve (see Fig.2 and Table 1). The 
agreement is fair at high magnetic fields, while at low magnetic fields is not good. We 
comment that the theory of the full GL model (higher Landau levels included) beyond 
Gaussian approximation is required at low magnetic fields. 



V. SUMMARY 

The problem of calculating the fluctuations effects in the framework of the Ginzburg - 
Landau approach is discussed and essentially solved in the LLL limit. The melting line 
location is determined and magnetization and specific heat jumps along it are calculated . 
The magnetization of liquid is larger than that of solid by 1.8% irrespective of the melting 
temperature, while the specific heat jump is about 6% and decreases slowly with 
temperature. The LLL results for melting line, magnetization jump, specific heat jump 
and the magnetization curves are in good agreement with experiments and Monte Carlo 
simulations. The leading corrections due to higher Landau levels have been also calculated 
and lead to agreement with experimental results in the whole region in which the GL 
approach is applicable. 
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Figure captions 

Fig. 1 

Free energy of solid (line) and liquid (dashed line) as function of the scaled temperature. 
The solid line ends at a point (dot) indicating the loss of metastability. Inset shows a tiny 
difference between liquid and solid near the melting point at a™ = —9.5. 

Fig. 2 

Comparison of the experimental melting line for fully oxidized YBa2Cu^0'j with our 
fitting. Inset on the right shows the relative universal magnetization jump of 1.8% (line) 
and experimental results for fully oxidized YBa2CuzO'j P (rhombs) and optimally doped 
untwinned YBa2Cu'iOT-s @] (stars). Inset on the left shows the relative nonuniversal 
specific heat jump (line) and experimental results for optimally doped untwinned 

Fig. 3 

Comparison of the theoretical magnetization curves (lines) of fully oxidized YBa2Cu^0-j 
utilizing parameters obtained by fitting the melting line on Fig. 2 with torque 
magnetometry experimental results (dots). Arrows indicate melting points while at low 
magnetic field the experimental data start from the point in which the magnetization is 
reversible indicating low disorder. 
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